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摘 要 : 在 正 交 闫 分 复 用 技术 和 双向 通信 下 ， 研 究 了 携 能 通信 网 络 中 的 资源 分 配 问 题 。 提 出 了 一 种 面向 子 载波 和 传 
输 功 率 的 联合 分 配 算法 ， 该 算法 不 依赖 于 现 有 的 功率 分 割 或 时 隙 切换 机 制 ， 可 以 在 不 同 的 子 载波 上 分 别传 输 信 息 和 
能 量 ， 简 化 了 携 能 机 制 设 计 。 此 外 ， 考 虑 到 实际 能 量 收 集 电 路 的 饱和 特性 ， 该 联合 算法 基于 非 线 性 能 量 采 集 模型 ， 
采用 拉 格 朗 日 对 偶 和 次 梯度 方法 ， 研 究 了 在 满足 用 己 传 输 速 率 阀 值 的 前 提 下 ， 使 系统 的 总 能 耗 最 小 化 的 问题 ， 降 低 
网 络 能 量 开销 。 数 值 仿真 实验 证 实 了 此 联合 分 配 算 法 的 有 效 性 和 能 量 高 效 性 。 
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Abstract: This paper investigated the resource allocation problem of the synchronous wireless information and power transfer 

(SWIPT) two-way communication networks based on orthogonal frequency division multiplexing. This paper proposes a joint 

subcarrier and power allocation (JSPA) algorithm, which is neither dependent upon the power splitting nor the time switching 

schemes and simplifies the design of SWIPT principles, transfers information and energy using different subcarriers. 

Considering the saturation characteristic of actual energy harvesting circuit, the joint design algorithm, which is developed by 

using Lagrange dual and subgradient methods and thus is based on the nonlinear energy harvesting model, is formulated to 

minimize the energy consumption of networks under condition that the required rate of the destination is satisfied. The 

numerical simulation results validate both the JSPA’s effectiveness in the joint allocation of network resources and its high 

energy efficiency. 
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评价 该 类 型 网 络 性 能 的 指标 之 一 。 在 满足 服务 质量 的 前 提 下 ， 


0 言 = bb E S7 Ap EN ALOA Z S 
JE 最 小 化 网 络 整 体能 量 开 销 是 提升 能 量 效率 的 有 效 途 径 。. 为 此 ， 
携 能 通信 (Simultaneous wireless information and power ”需要 从 算法 层面 着 手 ， 高 效 地 管理 与 分 配 网 络 可 用 资源 。 


transfer , SWIPT) 技 术 通 过 能 量 收集 (energy harvesting, EH) 电 
路 将 周围 环境 中 的 射频 (radio frequency, REF) 信号 转换 为 可 用 
的 直流 (direct current, DCO) 电源 趾 ， 为 能 量 受 限 的 传 感 占 、 移 
动 设备 、 低 功 耗 计算 机 等 终端 设备 供电 ， 在 传输 信息 的 同时 
解决 无线 通 信 网 络 中 设备 的 能 量 供应 问题 让。 据 测 算 ， 在 与 


Varshney! F 2008 年 首次 提出 了 信息 与 能 量 同时 传输 的 
策略 ， 访 策略 分 析 了 网 络 性 能 的 上 限 ， 但 由 于 射频 能 量 收集 
过 程 会 破坏 射频 信号 所 包含 的 数据 信息 ， 因 此 该 策略 未 能 实 
现 。 之 后 ，Zhou Xuniog 提 出 了 一 种 实际 的 信息 与 能 量 接收 机 
结构 ， 据 此 衍生 了 功率 分 割 (power splitting, PS) 和 时 际 切 换 


工作 于 GSM 频段 (1805~1880 MHz) 的 射频 源 距 离 50 米 处 的 
出 频 功率 密度 约 为 1mW/m?， 那么 以 100cm? 大 小 的 携 能 设备 
为 例 ， 它 将 最 多 收集 到 数 十 微 瓦 的 可 用 功率 ”。 这 对 于 超 低 
功率 需求 的 无 线 设 备 而 言 ， 已 具备 维持 其 间歇 应 用 的 基本 要 
求 。 所 以 ， 在 无 线 传 感 器 网 络 、 物 联网 等 5G 低 功 耗 受 限 网 
络 中 ，SWIPT 能 够 为 移动 用 户 提 供 了 极 大 的 便利 ， 已 引起 相 
天 学 者 的 关注 入。 在 无 线 携 能 通信 网 络 中, 能 量 效率 是 用 于 
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(time switching, TS) 两 种 通信 机 制 " 9， 为 携 能 通信 的 实际 应 
HAE S hiho 

基于 PS 通信 机 制 ， 文 献 [9] 研 究 了 OFDM 携 能 通信 中 继 
网 络 中 的 目 适应 传输 柑 式 转换 协议 。 该 协议 根据 源 剖 节操 与 
终 问 节点 不 同 的 通信 方式 ， 提 出 了 传输 模式 转换 策略 。 该 策 
略 使 用 拉 格 朗 日 对 偶 方 法 ， 求 解 具有 模式 转换 特点 的 非 串 资 
源 分 配 问题 。 基 于 TS 通信 机 制 ， 文 献 [10] 在 SWIPT NOMA 
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网 络 中 ， 以 最 大 化 能 量 效率 为 目标 ， 利 用 Dinkelbach 方法 ， 
提出 了 面 癌 传输 功率 和 TS 因子 的 双 层 资源 分 配 联合 优化 算法 。 
为 简化 接收 机 的 设计 复杂 度 ， 进 一 步 提 升 网 络 性 能 ， 一 
些 新 型 的 携 能 通信 策略 被 提出 ， 开 辟 出 信息 与 能 量 同 传 的 新 
途径 。 文献 [11] 提 出 了 新 型 的 接收 机 模式 切换 策略 , 其 中 多 天 
线 发 射 机 采用 传统 的 随机 波束 成 形 技术 发 送 公 共 信 息 。 每 个 
单 天 线 接收 机 在 未 知 信道 状态 信息 条 件 下 ， 根 据 无 线 信 道 的 
时 间 波 动 性 自主 切换 信道 ， 从 而 实现 回收 能 量 的 最 大 化 。 文 
献 [12] 针 对 SG 无 线 携 能 通信 网 络 研 究 能 量 效率 最 大 化 问题 ， 
并 提出 基于 OFDM 的 信息 与 能 量 传输 策略 , 对 网 络 资源 进行 
联合 优化 , 有 效 简 化 了 接收 机 的 设计 复杂 上 度 。 文 献 [13] 在 衰落 
信道 下 提出 了 动态 模式 切换 的 携 能 通信 方案 ， 以 最 大 化 用 户 
收集 能 量 为 目标 ， 上 自 适 应 调节 最 优 传输 功率 分 配方 案 。 

正 交 频 分 复 用 (orthogonal frequency division multiplexing 
OFDM) 技术 能 够 提高 网 络 频谱 效率 和 抗 多 径 干 扰 能 力 ， 有 
效 利 用 子 载波 的 传输 状态 ， 进 而 适 配 传输 功率 ， 能 够 显 音 降 
低 网 络 能 耗 , 提升 能 量 效 率 。 文 献 [14] 在 不 完美 信道 状态 信息 
KIFT, 研究 了 OFDM 携 能 通信 系统 中 的 资源 分 配 与 功率 控 
制 问题 。 文 献 [15] 在 下 行 正 交 频 分 多 址 网 络 中 , 为 满足 用 户 可 
乱 、 低 时 延 的 通信 需求 ， 提 出 了 一 种 基于 逐次 凸 优化 的 次 优 
资源 分 配 算法 ， 提 高 了 网 络 吞 吐 量 。 文 献 [16] 研 究 了 基于 
OFDM 技术 的 多 用 户 携 能 通信 系统 ， 访 系统 以 最 大 化 “和 速 
率 ” 为 目标 , 使 用 仿 禁 算法 实现 网 域 资 源 的 初始 分 配 。 虽然 它 
通过 任意 交换 两 个 用 户 的 子 载波 ， 提 高 了 中 晰 性 能 ， 但 随 兰 
子 载波 数量 的 增加 ， 算 法 复杂 上 度 也 将 大 大 增加 。 因 此 ， 与 传 
统 OFDM 通信 网 络 不 同 ， 基 于 OFDM 的 携 能 通信 网 络 ， 其 
性 能 将 同时 受到 用 户 数量 、 子 载波 峰值 功率 ， 以 及 信息 和 能 
量 传 输 功 京 水 平 的 影响 ”着 。 在 现 有 功率 分 配方 式 中 ,平均 功 
率 分 配 算 法 虽然 计算 复杂 上 度 低 ， 但 它 忽 视 了 信道 质量 ， 因 此 
当 子 载波 的 质量 较 差 时 ， 降 低 了 网 络 的 否 吐 量 。 

上 上述 相关 工作 大 多 是 基于 线性 EH 模型 来 分 析 网 络 性 能 
指标 。 然 而 近期 研 表明 清光 ， 基 于 二 极 管 的 EH 电路 的 RF- 
DC 转换 效率 是 非 线 性 的 , 而 不 是 线性 的 , 即 当 接 收 功 率 接近 
或 达到 饱和 国 值 时 ， 可 收集 的 射频 能 量 不 随 射 频 功 率 的 增 大 
而 增加 ， 有 共有 饱和 特性 。 因 此 ， 使 用 传统 的 线性 EH 模型 对 
系统 性 能 进行 分 析 和 设计 ， 可 能 会 导致 输出 不 准确 和 系统 性 
能 的 损失 。 文 献 [21] 提 出 的 符合 EH 电路 实际 工作 特性 的 非 
线性 EH 模型 ， 考 虑 了 射频 能 量 采 集 过 程 中 的 饱和 特性 ， 使 
得 网 络 性 能 的 评价 将 更 准确 , 共有 更 好 的 性 能 。 文 献 [22] 提 出 
了 一 种 简化 的 非 线 性 EH 模型 ， 在 可 操作 性 和 准确 性 之 间 取 
得 折 中 。 

本 文 所 研究 的 场景 为 基于 ODFM 技术 的 双 问 携 能 通信 
网 络 。 在 下 行 链 路 上 ,， 源 端 将 可 用 的 OFDM 子 和 载波 划分 为 两 
种 类 型 ， 分 别 用 于 射频 能 量 传输 和 信息 传输 。 因 此 ， 终 端 接 
收 机 无 须 额 外 设计 功率 分 割 或 时 隐 切 换 机 制 ， 简 化 了 携 能 机 
制 设计 。 由 于 射频 EH 电路 存在 饱和 特性 , 基于 非 线性 EH 模 
型 , 本 文 提出 了 子 载波 和 传输 功率 的 联合 分 配 Qoint subcarrier 
and power allocation, JSPA) 算 法 。 数 值 实验 结果 证 实 ， 在 满足 
用 户 数据 速率 需求 的 前 提 下 , JSPA 算法 能 够 得 出 子 载波 集合 
的 最 优 划 分 方案 ， 揭 示 子 载波 与 功 座 、 传 输 速 率 和 能 量 回收 
之 加 的 内 在 联系 ， 有 效 减少 网 络 能 量 开销 。 


1 ”网 络 模型 和 问题 搬 述 


1.1 网 络 模型 

考虑 如 图 1 所 示 的 OFDM 携 能 双 同 通信 网 络 , 该 网 络 由 
具备 单 天 线 的 源 问 S 和 具有 能 量 收集 能 力 的 终端 DHAR. 
端 贡 点 的 全 部 能 量 来 源 于 收集 源 端 发 射 的 RF 信号 ， 其 将 接 
收 到 的 射频 信号 用 于 能 量 收集 和 信息 解码 ， 并 使 用 回收 的 射 


Z, F: 双向 携 能 通信 网 络 中 子 载波 和 功率 联合 分 配 
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一 > 源 端 向 终端 传输 信息 的 子 载波 
-一 -> 源 端 向 终端 传输 能 量 的 子 载波 


一 一 终端 同 源 端 通信 时 际 ”所 一 一 终端 向 源 端 传输 信息 的 子 载波 
T, 


T y Al Al bi > Zd D 
一 全 一 源 端 向 终端 通信 时 阶 


图 1 OFDM 双向 携 能 通信 系统 模型 

Fig. 1 OFDM SWIPT based tow-way communication system model 

FE Wa Sing BY ZA ig AY FE Ss STE, WEIN RA, AS 
宽 集合 被 均等 划分 为 N 个 正 交 子 载波 N= 全 .…,n.…N} ,其 
中 ， 第 7 个 子 载波 的 信道 增 蔡 为 h ， 源 病 为 子 载波 4 分 配 的 
BN TED p, 。 在 终端 到 源 端 的 单 跳 阶 段 ， 通 信 时 际 为 区 ， 
可 用 市 宽 集 合 MY 被 均等 划分 为 AL= 全 …m.…M}， 其 中 ， 8n 
为 第 mw 个 子 载波 的 信道 增益 ，4, 为 终端 发 射 功率 。 中 各 
分 别 表示 源 端 与 终端 均 服 从 正 态 分 布 的 加 性 高 斯 白 噪 声 。 源 
端 与 终端 均 采 用 基于 OFDMA 接收 机 的 带 通 滤波 器 ， 并 已 知 
信道 状态 信息 ， 可 以 实现 在 不 同 的 子 载波 上 分 别传 输 信息 和 
功率 请]。 此 外 ， 由 混 频 器 、 滤 波 器 和 数 模 转换 器 等 模块 消耗 
的 电路 功率 与 发 射 功率 相 比 非常 微小 ， 本 文 忽略 不 计 。 
1.2 问题 描述 

源 端 S 通过 OFDM 子 载波 集 N 传输 信息 和 能 量 到 终端 
D, REES N 补 划 分 为 两 个 子 集 ,好 GG 和 Cs(G,CcN,， 
GCIG=N , GAG=D), 如 图 1 所 示 , FRG 中 的 子 载波 用 
于 信息 传输 , 其 子 载波 i 上 的 传输 功率 记 为 vi (EG) FÆ C 
中 的 子 载波 用 于 能 量 收集 ， 其 子 载波 /上 的 传输 功率 记 为 对 
(jeCs) 下 行 通信 和 链 路 上 可 实现 的 数据 传输 速率 为 


pih, 
及 。 = Í [p= 
co Dilog| 1+ 2) (1) 
源 端 因 传 输 能 量 和 信息 到 终端 而 产生 的 能 耗 为 
N 
Ws=T YB (2) 
n=l 
终端 采集 到 的 能 量 表 示 为 


下 = 站 ACen 
之 /7 (Fin) (3) 


其 中 ， fi() 为 文献 非 线 性 能 量 转换 模型 ， P= ph top. Je 。 
在 能 量 采 集 电路 的 输入 射频 功率 范围 内 ， 能 量 转换 效率 
通 第 先 随 输 入 功率 的 增 大 而 增加 ， 到 达 国 值 后 的 转换 效率 随 
着 输入 功率 的 增加 而 下 降 。 根 据 能 量 采 集 过 程 中 转换 效率 的 
非 线 性 ， 能 量 采 和 集 模型 可 近似 为 ?231;: 
fu (Pa) =@,P2 + oP, + 0 (4) 
其 中 ，% 、2% Mla, 为 与 能 量 采 集 过 程 有 关 的 参数 ， 由 终端 
天 线 和 能 量 采 集 电 路 的 设计 而 定 。 
在 终端 到 源 端 的 上 行 链 路 上 ，OFDM 子 载波 均 用 于 信息 
传输 ， 上 行 链 路 的 可 达 速 率 Ros 为 


Reyye 2 log( 1+ 28e (5) 
me At Nn 
24 im Al FE Hn es JS BN im TT P= AE A BE EA 
Wp =T,’ „an (6) 


EZA RENAR, Mim Ee ae ERE S 
运行 。 以 最 小 化 源 端 能 耗 为 目标 ,采用 非 线 性 能 量 采集 模 型 ， 
将 源 端 和 终端 的 速率 需求 和 能 量 收集 门限 纳入 约束 ， 联 合 优 
化 子 载波 和 传输 功率 ， 得 到 如 下 优化 问题 PA: 
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min TYP 7 fie Ee gs, N 11 
{Pn,dm,C7 Ce} i i ( a) N, + q, 8 m í ( ) 
s.t. Soe 1422) >0, (7b) af Ry -Slog(1+ £28 ]) 0 (12) 
icC oF m= Mn 
T, >) fu (ph; +0})-Wp 20, (7c) ada ie (13) 
i 得 到 In 解 的 表达 式 为 
Silos 1+ 48 |-R, 20, (7d) an (2-2) me M (14) 
m=1 m Em 
Pin 2 Pror Pin S Pos (7e) 其 中 ， xt =max(0,x) 。 
P, 20,.ne MV (7f) 将 式 (14) 代 入 式 (10)， 重 写 子 优化 问题 P2 的 拉 格 明日 函数 为 
Gn 2 0,m E€ AL (7g) 


其 中 , 约束 条 件 (7b) 和 (7d) 分 别 为 源 端 和 终端 用 户 的 最 小 数据 
传输 速率 办 什 要求， 约束 条 件 (7c) 表 示 源 端 发 射 的 传输 功率 
应 满足 能 量 收集 门限 , 即 提供 给 终端 进行 数据 传输 所 需 的 能 量 。 
约束 条 件 (7e) 表 示 能 量 收集 电路 可 接受 的 的 输入 功率 范围 , 约束 
条 件 (7D 和 (7g) 要 求 源 端 和 终端 的 发 射 功率 为 非 负 值 。 


2 子 载 波 和 功率 联合 分 配 算 法 


上 述 优化 问题 P1 中 ,目标 函 数 (7a) 和 约束 条 件 (7b)~(7d) 
均 为 凸 函 数 ，(7e)~(78) 为 线性 约束 条 件 ， 且 存在 多 个 待 确定 
的 集合 变量 : P= pPlieC.ieG} , O={q,|meAQ , 
el | GEN e rG O e A 

若 利 用 穷 举 搜索 法 根据 优化 问题 Pl 寻找 子 载波 的 最 优 
分 配方 案 ， 其 算法 复杂 上 度 会 随 独 子 载波 数量 的 增加 呈 几 何 级 
数 式 大 大 增加 。 通 过 观察 双 回 通信 网 络 的 特点 和 问题 结构 ， 
原 优化 问题 Pl 中 与 优化 变量 8 相关 的 约束 条 件 只 有 (7d) 和 
(7g), 而 且 该 约束 条 件 决 定 了 终端 进行 数据 传输 所 需 的 功率 ， 
与 优化 问题 中 的 其 他 约束 条 件 和 目标 函数 无 耦合 。 因 此 ， 为 
了 降低 算法 复杂 上 度 ， 原 优化 问题 Pl 可 分 解 为 两 个 独立 的 子 
优化 问题 P2 和 P3。 子 优化 问题 P2 以 最 小 化 终端 能 耗 为 目 
标 ， 沿 用 原 问 题 中 的 (7d) 和 (7g) 作 为 约束 条 件 ， 其 表达 式 为 


M 
minT, $ q, 
Im m=l 


M 
S. t. > os| + dnn ) — Rp 20, (8) 
m=] Nn 


I 


n 2 0, m E€ M. 


该 双向 携 能 通信 网 络 中 ， 终 端 可 用 的 能 量 来 源 于 从 源 端 
收集 到 的 射频 能 量 。 子 优化 问题 P2 以 最 小 化 终 站 能 耗 为 目 
标 ， 当 终端 能 耗 最 小 时 ， 源 端 能 耗 也 降 至 最 低 。 因 此 ， 求 解 
P2 问题 所 得 到 的 最 优 解 集 ， 是 原 问题 Pl 解 集 的 子 集 。P1 经 
剥离 子 问 题 P2 后 得 到 的 新 问题 ， 可 记 为 优化 问题 P3: 


min T, > Pn 
n=l 


{PC} 


Th. 
S.t. Yitoo(1+ 24 )- r > 0, 
icC, Oi 


T, >, fr(pFh;+o?)—Wo >0, 
JECF 


(9) 


pi 20,1€C, 
Pin Z Piws Pin S Pup 
2.1 网 络 模型 优化 终端 功率 分 配 
针对 子 优化 问题 P2， 设 定时 际 ， 为 使 终 背 能 量 消耗 最 
小 化 ， 需 要 求解 终 问 的 最 优 功率 分 配方 案 。 为 此 ， 构 建 子 问 
题 P2 的 拉 格 并 日 函数 : 


M M 
LC(qn, Kn, 1) 一 Gs ~ > eg: + 
m=l m=l 


af Re -$ tog{ 14 dngn ) (10) 


2 
Nin 


其 中 ， kna (K, 20, me AD) 和 4 (4>0) 为 拉 格 朗 日 乘 子 。 
结合 子 问 题 P2 的 约束 ，Karush-Kuhn-Tucker 条 件 可 表示 为 


M g M n2 
L(A)=A(M kr as a (15) 
m=1 m m=1 Sm 


此 时 ， 式 (15) 为 关于 2 的 四 函数 ,， 令 其 一 阶 导 数 为 零 ， 得 
到 对 偶 变 量 4 的 最 优 值 为 


M 
jean ae (16) 


将 最 优 乘 子 生 代入 去 (14)， 可 得 终端 友 大 功率 的 最 优 值 为 


M 
CRp-2 M n? 
m=| m — 


Im =(e +, m E€ M (17) 
Em 

子 优化 问题 P2 eA A im BE FE Wo =LA) ce VG m AA 
能 量 收集 门限 ， 代 入 子 问题 P3 的 能 量 采 集约 束 条 件 


TY) fu (pfl +07) -Wo 20 后 更 新 子 优化 问题 P3。 


2.2 优化 拉 格 明日 对 偶 变 量 

通过 求解 子 优化 问题 P2, 得 到 终端 能 耗 最 小 化 的 最 优 值 
Wb。。 子 优化 问题 P3 以 最 小 化 源 端 能 耗 为 目标 ， 以 源 端 最 小 
数据 传输 速率 需求 和 最 小 能 量 采 集 门 限 为 约束 条 件 ， 设 定时 
阶 五 =1。 利 用 拉 格 天 日 乘 子 法 ， 得 到 子 优化 问题 P3 的 拉 格 
明日 函数 如 下 : 


N pih, 
LP, C)=}_ p, mG -os 路 
n=l ieC7 i 


(18) 

pz (w 一 2, (œ (ph yh? +æ přh, + Qs ) 

其 中 ， 2 和 有 (6>0， 有 >=0) 为 拉 格 朗 日 乘 子 。 

子 优化 问题 P3 的 拉 格 明日 对 偶 函 数 为 
g(B)= min L(P,C) (19) 
EF, B=(8.4) 。 得 到 子 问 题 P3 的 对 偶 优 化 问题 为 

max g(B) m 
s.t. B20 ou 


fE H Ui tks PSE SEAS ORS APPS B Ted LY ee B,C PE 


I 


AB, =Wp — > (a (pF)? h? + a, pFh; +a) 
JeCE 


KHR, =b), BIRI Bo =P + db 逐渐 更 新 。 根据 减 小 
PEDT, AREE v =w/t 。 其 中 ，。 为 常数 ， 次 梯度 方法 
可 保证 对 偶 变 量 收敛 于 最 优 解 广 =4p) PA, 
2.3 源 病 功率 分 配 
根据 对 偶 函 数 式 (19)， 在 确定 的 对 偶 变 量 =A) F, 
求解 源 端 发 射 功 率 分 配 最 优 解 p 。 对 于 给 定 的 子 载 
WL AG C={G,G) ， 可 对 式 (18) 所 示 的 拉 格 并 日 函数 其 优化 变 
量 p, 求 偏 导 ， 得 (22) 式 : 
OLC, pi) 1 _ Ph 
Op} o? + pih, 


OL(C, p7) 
= a =1— 201p,p?h? 一 fh, 


AB, = Rs -Yitog{ 1+ 2) 
sa 5 (21) 


(22) 


J 


根据 Karush-Kuhn-Tucker 条 件 ， 拉 格 表 日 函数 的 偏 导数 
在 最 优 解 处 为 零 。 因 此 ， 源 端 功率 最 优 解 Z 为 


* o? j 
p= [2 -至 (23) 
h, 
Prins PFA; +07 < Prow 
1—a,f,h, i 
E j 
a | 20 Bh; 24) 


Dar Di Nye tO; > Pu 
其 中 ， Prin =0， Pam = (pw -07)/h 
2.4 ”优化 子 载波 分 配 
根据 对 偶 函 数 式 (19)， 在 确定 的 对 偶 变 量 B=(4B) F, 
将 求 得 的 源 端 最 优 发 射 功 率 Z = (7, Pe} 代入 (18) 式 所 示 的 拉 
格 表 日 函数 ,对 其 整理 并 化 人 简 为 (25)(26) 式 , 求解 子 载波 分 配 
最 优 解 C GG} 。 


ih, 
L(C) = BR; + Bar +3 or 一 局 oe0+ 205) 
iel; i 
> (pF -p (a (ph, +07) +a, (pF hj +0})+a;))= 


jel, 
BRs + >》 (vp -Ba (pPh, +07) +a, (pFh +02)+a)- 
JECF 


(25) 
mh N r 
pj + logd+ L n+ Bar +> ve ~ B log + a | = 
DG + AR, + Br +S ve - p, logd + 2% | 
jelz n=1 O; 
KE, B=) 。 得 到 子 问题 P3 XBA A 
G; = př —a,B, (ph, +03) - 
(26) 


i : pih; 
Anh (PDF h, +07)—Qp, =p + B, Og og? 
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eH, op-[6-Z] , jea. 同时， 为 简化 计算 ， 令 
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在 式 (25) 中 ,只 有 第 一 项 受 载波 集合 A 划分 结果 的 影 
Mel, 其 他 项 均 为 单数 。 因 此 , 根据 仪 包 含 优化 变量 G 的 式 (26)， 
可 得 到 集合 G 的 最 优 解 Ce ， 


C: =argmin > Gř. 
E 8 区 2 j (27) 


在 载波 集合 V 的 划分 结果 ( 即 G 和 Gi ) 未 知 条 件 下 ， 为 
使 子 问题 P3 的 拉 格 朗 日 对 俩 函数 最 小 , 应当 将 子 载波 j 在 集 
BM PR, BID 28 ERE Pr A Pr RAS (26). AIN 
(26) 的 函数 值 为 负 , 则 将 该 子 载波 7 划分 到 集合 GP RA, 
所 有 能 使 26) 式 函数 值 为 负 的 子 载波 所 构成 的 集合 ， 即 为 
用 于 能 量 传输 子 载波 的 最 优 解 Ceo IR UG HN, 可 得 到 
用 于 信息 传输 的 最 优 载 波 集合 0 为 

G =G. (28) 

初始 化 对 偶 变 量 BO), ERR K=O, MERTER A 
功率 联合 分 配 的 最 优 解 2 = PP AR CIE) o RE A All 
”代入 式 (21)， 可 得 P3 对 偶 优 化 问题 的 次 梯度 Ag ， 更 新 
k=k+1 ， 进 而 根据 次 提 梯 度 下 降 方 法 得 到 BK). 
B-B- RENAR, 根据 BC) 重新 计算 A AC ， 
否则 ,算法 结束 , 输出 P3 目标 函数 的 最 优 值 。 初 始 优化 问题 
Pl 的 求解 方法 可 归纳 为 联合 子 载波 和 功率 分 配 Joint 
subcarrier and power allocation, JSPA) 算 法 ,算法 流程 图 如 图 2 
所 示 。 

算法 1 联合 子 载波 和 功率 分 配 算法 

a) 设置 和 迭代 次 数 上 =0， 初 始 化 及 ,RdTDPMN ; 

b) 根据 式 (10) 计 算 和 ， 根 据 式 (17) 计 算 (an Wo) 

c) 初始 化 4.0) =1,4,0)=1 
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d) 根据 式 (23) 和 (24)， 计 算 功 率 分 配 最 优 解 p 和 Pp 

e) 根据 式 (27) 和 (28)， 得 到 子 载波 分 配 最 优 解 必 和 

f) 更 新 k=k+1 ， 根 据 式 (21) 更 新 对 侦 变 量 A0,A0) 

g) WMR|BO-Bk-D)kKelBwO-Bk-D\k<e, BRIE, 输 
HE A, LN ALG GPP PE hd; AR, BRER 四 继续 
循环 。 

提出 的 JSPA 算法 输入 量 包 括 源 端 和 终端 的 最 小 数据 传 
fen OR SRA Rs AU R ， 源 端 与 终端 之 间 的 传输 距离 4 ， 上 行 和 
下 行 链 路 划分 的 子 载波 数量 M 和 N ,以 及 上 下 行 链 路 的 传输 
IN BT 和 。 算 法 输出 结果 包括 下 行 链 路 分 别 用 于 信息 传输 
和 能 量 传输 的 子 载波 分 配方 式 G 和 把 ， 以 及 上 下 行 链 路 各 
子 载波 上 的 功率 分 配方 案 grmesM M pi,neN 。 

优化 问题 P1 通过 分 解 等 价 转换 为 P2 和 P3 两 个 子 优化 
问题 问题 。 其 中 , 子 优化 问题 P2 为 凸 优化 问题 ， 通 过 拉 格 歧 
日 乘 子 法 可 得 到 其 最 优 解 ， 而 子 优 化 问题 P3 约束 条 件 中 存 
在 优化 变量 的 相互 耦合 ， 为 非 凸 优化 问题 ， 为 了 得 到 其 最 优 
解 ， 需 要 先 将 P3 转换 为 上 四 优化 问题 。 根 据 目 优化 理论 , 无 论 
原 优化 问题 是 否 为 凸 优 化 问题 ， 其 拉 格 衣 日 对 侦 问 题 一 定 为 
凸 优化 问题 P3， 将 非 凸 的 P3 优化 问题 转换 为 凸 的 对 偶 优 化 
问题 式 (20)， 式 (20) 的 优化 变量 满足 time-sharing 条 件 P4， 则 
式 (20) 与 P3 优化 问题 之 间 为 强 对 偶 关 系 ， 因 此 式 (20) 的 最 优 
解 就 是 P3 的 最 优 解 ，P2 和 P3 的 最 优 解 保证 了 提出 的 JSPA 
算法 求 得 的 最 优 解 为 原 优化 问题 P1 的 最 优 解 。 


1 | 初始 化 RoR; 、M、N、d 等 相关 参数 


采用 拉 格 朗 日 乘 子 法 ， 解 决 子 
优化 问题 P2， 利 用 KKT 条 件 和 


点 求 时 公式 求 得 对 偶 变 量 最 优 解 X 
, i 
”根据 式 (7 计算 终端 传输 功率 
最 优 值 4;, 和 终端 最 小 能 耗 W 
a ae oe a + a a ee A 
CE A 


针对 子 优化 问题 P3， 采 用 拉 格 朗 日 
对 侦 方 法 ， 得 到 其 对 侦 优 化 问题 


} 
利用 KKT 条 件 ， 在 确定 的 对 偶 变量 6 
下 ， 根 据 式 (23) 和 (24) 计 算 功率 分 配 
最 优 值 p; A p7 


在 已 知 的 对 侦 变 量 和 优化 变量 {p;, p;7 
ge | 下 ， 采 用 换 元 法 重 写 子 优化 问题 P3 的 

5 | 拉 格 朗 日 函数 ， 并 根据 式 (27) 和 (28) 计 
算 子 载波 分 配 最 优 值 GM CG 

采用 次 梯度 方法 根据 式 

(21) 更 新 对 偶 变 量 {5., 9, } 


> 


不 
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TE 


HFA fil EPS BOIL 


2 子 载波 和 功率 联合 分 配 算法 流程 图 
Fig.2 The flow chart of joint subcarrier and power allocation algorithm 
本 算法 在 拉 格 明日 对 偶 变 量 求解 步骤 中 采用 的 次 梯度 算 
法 ， 其 算法 复杂 度 为 CU PI, AEL 为 对 侦 变 量 个 数 ，x* 为 
非 负 整数 。 由 于 本 算法 只 涉 hb 和 两 个 对 侦 变 量 ， 因 此 ,该 
步 又 的 算法 复杂 度 为 0(2*)=00) 。 在 子 载波 分 配 求解 步骤 中 
算法 复杂 度 和 子 载波 数量 N 成 正比 ， 该 步骤 的 算法 复杂 度 为 

O(N) ， 因 此 算法 整体 计算 复杂 上 度 为 O04dxN)=0(N)。 
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本 节 分 析 和 验证 子 载波 和 传输 功率 联合 分 配 算法 的 有 效 
性 。 在 微波 功率 传输 可 行 环境 中 ， 源 端 发 射 功率 与 终端 接收 
功率 间 的 关系 满足 式 (29)1， 

P. _ AA, 


P, A2d? 
其 中 ，4 是 波长 ，4 是 源 端 到 终端 的 传输 距离 ，4 和 4 分别 
是 发 射 天 线 阵 列 和 接收 天 线 阵 列 的 总 孔径 ，2Z 为 复 高 斯 随机 
变量 ， Z ~ CN (1,0.2) .本 文中 Pe i LAF = 2.4GHz , v, =3x10° m/s ， 
d=2m, A,=0.1m?, A, =0.05m 281， 子 载波 数量 W=M=4， 满 足 
源 端 和 终端 服务 质量 需求 的 数据 传输 速率 为 
Rs =R, =1bit/sec/ Hz ， 源 闫 用 于 信息 传输 信道 上 的 噪声 
o?7=0.001(ieG)， 用 于 能 量 传 输 信 道上 的 噪声 07=0( jeG)， 
终端 受到 的 噪声 干扰 na =0.001( me A4 )。 能 量 收集 电路 所 允许 
的 输入 功率 范围 为 B=103WW , Pp =32004W [2], 

图 3 给 出 了 在 线性 EH 模型 和 非 线 性 EH 模型 下 ， 源 端 
最 小 能 耗 与 能 量 收集 门限 的 关系 。 当 线性 能 量 收集 (Linear EH, 
LEH) 模 型 运用 于 子 载波 和 功率 联合 分 配 算法 时 , 终端 将 根据 
线性 能 量 转换 效率 获得 可 用 能 量 。 由 图 3 可 知 , “4 LEH 模型 
中 转换 效率 8 为 0.2 和 0.4 时 , 源 端 最 小 能 耗 高 于 使 用 非 线性 
能 量 收集 (Non-Linear EH, NLEH) 模 型 时 的 最 小 能 耗 。 在 LEH 
算法 中 ， 念 真 结果 显示 终端 实际 采集 到 的 能 量 高 于 其 能 量 收 
集 门限 ， 说 明 功 率 与 子 载波 在 分 配 机 制 设计 上 不 合理 ， 加 剧 
了 源 端 能 量 消耗 。 当 LEH 模型 中 转换 效率 8 为 0.6 时 ， 源 端 
最 小 能 耗 值 低 于 采用 NLEH 模型 时 的 源 端 能 耗 。 这 表明 LEH 
模型 设置 过 高 的 能 量 转换 效率 值 了 会 导致 分 析 结 果 过 于 理想 
偏离 了 实际 ， 如 果 采 用 该 方案 分 配 网 络 资源 ， 网 络 能 耗 虽 少 
但 由 于 未 考虑 EH 电路 的 饱和 特性 ， 其 实际 收集 能 量 并 达 不 
到 用 户 能 量 收集 门限 。 绽 上 所 述 , 采用 LEH 模型 分 析 携 能 通 
信和 网 络 性 能 时 ， 无 论 转 换 效率 8 设置 为 何 值 ， 其 分 析 结 果 终 
会 与 实际 网 络 性 能 有 所 信和 关 ， 而 基于 NLEH 模型 的 子 载波 和 
功率 联合 分 配 算法 能 够 根据 实际 EH 电路 的 非 线 性 特性 ， 合 
理 有 效 地 分 配子 载波 和 传输 功率 等 网 络 资源 ， 满 足 双 回 携 能 
网 络 对 数据 传输 速率 的 需求 ， 分 析 结 果 更 符合 实际 。 


12 (29) 


A LEH 模型 转换 效率 5=0.2 
0:035 f | 一 se 一 LEH 模型 转换 效率 6=0.4 

一 x 一 LEH 模型 转换 效率 =0.6 
| | 一 e 一 NLEH 模型 


源 端 最 小 能 耗 W _ (W) 


0 100 200 300 400 500 600 700 800 900 
终端 能 量 收集 门限 W (uW) 


3 源 端 最 小 能 耗 随 能 量 收集 门限 的 变化 


Fig.3 Minimal energy consumption of the source vs. threshold of 


energy harvested at destination 

4 显示 了 子 载 站 和 功率 联合 分 配 (joint subcarrier and 
power allocation, JSPA) 算 法 与 基准 算法 的 性 能 对 比 。 基 准 算 
法 包括 : 分 别 基 于 PS. TS 机 制 的 功率 分 配 算法 、 固 定子 载 
波 分 配 (fixed subcarrier allocation, FSA) 和 固定 功率 分 割 (fixed 
power splitting, FPS)。 在 基于 PS 机 制 的 功率 分 配 算法 中 ， 系 
统 接收 器 采用 功率 分 割 结 构 设 计 ， 各 个 子 载波 上 的 传输 功率 
通过 PS 因子 * 将 分 为 两 个 部 分 , x 部 分 功率 用 于 传输 信息 ， 
l-r 部 分 功率 用 于 传输 能 量 。 在 基于 TS 机 制 的 功率 分 配 算 
法 中 ， 系 统 接 收 顷 采用 时 隐 切 换 结 构 设 计 ， 终 姗 将 信号 传输 
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的 单位 时 隙 通过 TS 因子 吧 分 为 两 个 阶段 : 时 际 的 第 一 个 阶 
段 wz 用 于 能 量 收集 ， 时 际 第 二 个 阶段 4-o) 二 用 于 信息 解码 。 
在 FSA 算法 的 下 行 链 路 上 ， 两 个 子 载波 用 于 信息 传输 ， 另 两 
子 载波 用 于 射频 能 量 传 输 ， 而 源 端 与 终端 的 传输 功率 则 通过 
注水 法 加 以 分 配 。 在 FPS 算法 中 ， 源 端 与 终端 分 别 根据 子 载 
波 的 质量 为 其 分 配 传输 功率 ， 以 使 得 每 个 子 载波 上 的 数据 速 
RETEK BE. 

由 图 4 可 知 ， 本 文 所 提出 的 JSPA 联合 优化 算法 ， 在 源 
iy B/D BERET A, 整体 优 于 其 他 对 比 算法 。 痛 先 , 基于 PS 和 
TS 机 制 的 功率 分 配 算法 未 考虑 网 络 信道 的 可 分 配 性 , 因此 不 
能 根据 子 信 道 间 的 差异 特性 灵活 分 配 传输 功率 。 随 着 源 端 数 
据 速 率 靖 值 的 增加 ， 相 比 其 他 考虑 信道 分 配 的 算法 ， 该 两 算 
法 的 功率 分 配方 案 大 大 增加 了 源 端 的 能 耗 。 其 次 ，FSA 算法 
忽略 了 信息 传输 与 能 量 传输 的 不 同 特 性 ， 其 固定 不 变 的 信道 
划分 和 使 用 方式 ， 浪 费 了 源 问 能 量 。 然 后 ，FPS 算法 根据 速 
率 冰 值 为 子 载波 分 配 功率 ， 这 种 方式 无 法 实现 功率 与 子 载波 
质量 间 的 最 佳 匹 配 ， 过 度 消耗 了 源 端 能 量 。 可 预见 的 是 ， 当 
选用 频率 选择 性 衰落 信道 时 ， 子 载波 质量 的 差异 性 将 明显 增 
强 ，FPS 算法 中 的 源 端 能 耗 也 将 增 大 。 最 后 ，JSPA 算法 中 ， 
子 载波 按照 功能 (信息 传输 、 能 量 传输 ) 实 现 了 最 优 角色 划分 ， 
而 且 每 条 载波 上 的 传输 功率 值 也 最 佳 匹 配 于 载波 质量 和 载波 
HE. KE, JSPA 算法 能 够 在 满足 数据 速率 半 值 的 前 提 下 ， 
降低 源 端 能 量 开 销 。 
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4 源 端 最 小 能 耗 随 源 端 数据 速率 闪 值 的 变化 
Fig.4 Minimal energy consumption of source vs. threshold of 
required rate at destination 

子 载 波 和 功率 联合 分 配 算法 的 收敛 行为 如 图 5 和 6 所 示 。 
可 以 看 出 当 距 离 4 为 2m， 能 量 收集 门限 W 为 0.9mw， 最 小 
BIE Fan WRK BY AL Rs = [bit / sec/ Az WY, GEIR AL t 的 增 
加 ， 图 5 PAAR BAA He A 收敛 到 最 优 解 0.0025. 2 Wc st Hi 
最 优 解 18.15, 图 6 得 到 最 小 化 源 端 能 耗 优化 问题 的 最 优 值 。 
联合 观察 图 $ #6, TEER BON 20 左右 ， 源 端 总 能 耗 数 值 
基本 收 钱 ,体现 了 提出 的 子 载波 和 功率 联合 分 配 算法 的 有 效 
的 收敛 性 。 
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Fig. 5 Convergence behavior of Lagrange multipliers 4, and /, 
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a 的 传输 功率 不 能 达成 能 耗 最 小 化 的 目标 时 , JSPA 算法 将 动态 
一 # 一 子 载波 和 功率 联合 分 配 自 法 的 收敛 行为 调整 子 载波 角色 ， 划 分 更 多 的 子 载波 用 于 数据 传输 。 

在 源 问 能 量 收集 门限 (终端 最 小 能 耗 ) Wo 为 500W 、 
1500uW 和 2500uW 时 ，JSPA 联合 优化 算法 得 到 的 资源 分 配方 
案 如 图 10 所 示 。 由 图 10 可 知 ， 当 固定 传输 距离 和 数据 速率 
闵 值 时 ， 随 着 能 量 收集 门限 的 增 大 ， 为 维持 终端 增 大 的 能 量 
需求 ， 源 端 需要 分 配 更 多 数量 的 子 载波 用 于 射频 能 量 传输 ， 
导致 用 于 传输 数据 的 子 载波 的 数目 减少 。 综 合 观察 图 9 和 图 
10, 源 问 为 子 载波 分 配 的 功率 主要 用 于 能 量 收集 .这 是 因为 : 
a) 射 频 能 量 的 转换 效率 较 低 ， 需 要 更 多 的 传输 功率 才能 达到 


250 r 


BT aptamer o S a REWER I ER OE RAGER KED, DM FEZ E 
6 子 载波 和 功率 联合 分 配 算法 收敛 行为 率 即 可 满足 速率 需求 。 
Fig.6 Convergence behavior of joint subcarrier and power allocation algorithm 0.018 : 
终端 最 小 能 耗 ( 源 端 能 量 收集 门限 ) 与 数据 速率 阔 值 和 伟 soro | | mw’ 
输 距 离 4 的 关系 如 图 7 所 示 。 由 图 7 可 知 ， 随 着 数据 速率 阐 TH Fs 
值 的 增 大 ， 终 端 最 小 能 耗 随 之 增加 。 妆 传输 距离 固定 时 ， 速 ootz| 
率 阀 值 增 大 迫使 终端 消耗 更 多 功率 提高 可 达 数 据 速率 。 当 回 = oor} 
ERK BREET , 传输 距离 的 增加 降低 了 信道 传输 质量 , 此 时 ， È oos } 
2A inn ma 22 1B FE E L HS RE EE Y AE A r E FF AE FR E ER o ~ 0.008 | 
源 端 最 小 能 耗 与 能 量 收集 门限 (终端 最 小 能 耗 ) 和 传输 距 0.004 | AA 
离 4 之 间 的 关系 如 图 8 所 示 。 由 图 8 可 知 ， 当 数据 传输 速率 oooz| i, 
闵 值 和 能 量 收集 门限 固定 时 ， 随 着 传输 距离 的 增 大 ， 链 路 质 M BAY 
量 下 降 , 源 端 发 射 更 多 的 功率 用 于 信息 传输 导致 总 能 耗 增加 。 人 
当 数 据 传 输 速率 国 值 和 传输 距离 不 变 时 , 能 量 收集 门限 越 大 ， 图 9 子 载波 和 功率 分 配 随 源 端 数据 速率 闷 值 的 变化 
源 问 需要 消耗 更 多 能 量 用 于 能 量 传 输 ， 这 导致 忌 能 耗 增 大 。 Fig.9 Allocation of subcarriers and power vs. 
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Fig.7 Minimal energy consumption of the destination vs. threshold of Fig. 10 Allocation of subcarriers and power vs. 
required rate at destination and transmission distance threshold of energy harvested at destination 
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S| |- woow 功率 的 联合 分 配 算法 。 源 端 将 可 用 的 OFDM 子 载波 划分 为 两 
a 组 , 分 别 用 于 信息 传输 和 射频 能 量 传输 JSPA 算法 以 源 端 能 
= a, 量 开销 最 小 为 目标 ， 根 据 能 量 采集 电路 的 饱和 特性 ， 采 用 非 
要 线性 能 量 采 集 模型 ， 在 满足 速率 阔 值 要 求 的 前 提 下 ， 得 出 子 


载波 最 优 划 分 和 功率 匹配 结果 。 与 传统 携 能 通信 机 制 不 同 的 
是 无 须 功 率 分 流 器 和 时 际 切换 器 , JSPA 简化 了 携 能 机 制 设计 。 
TE a a oer ae 下 一 阶段 工作 欲 将 人 工 噪 声 引 入 多 输入 多 输出 的 双 工 
传输 距离 d (m) SWIPT 网络, 着 重 提升 用 户 信息 安全 性 , 最 大 化 网 络 能 量 效率 。 
8 源 端 最 小 能 耗 随 能 量 收集 门限 和 传输 距离 的 变化 参考 文献 ， 
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